Deciphering the Multi-state Conformational Equilibrium of HDM2 in the Regulation of p53 Binding: Perspectives from Molecular Dynamics Simulation and NMR Analysis.

HDM2 negatively regulates the activity of the tumor suppressor p53. Previous NMR studies have shown that apo-HDM2 interconverts between an "open" state in which the N-terminal "lid" is disordered and a "closed" state in which the lid covers the p53-binding site in the core region. Molecular dynamics (MD) simulation studies have been performed to elucidate the conformational dynamics of HDM2, but the direct relevance of the experimental and computational analyses is unclear. In addition, how the phosphorylation of S17 in the lid contributes to the inhibition of p53 binding remains controversial. Here, we used both NMR and MD simulations to investigate the conformational dynamics of apo-HDM2. The NMR analysis revealed that apo-HDM2 exists in a fast-exchanging equilibrium within two closed states, closed 1 and closed 2, in addition to a previously demonstrated slow-exchanging "open-closed" equilibrium. MD simulations visualized two characteristic closed states, where the spatial orientation of the key residues corresponds well to the chemical shift changes of the NMR spectra. Furthermore, the phosphorylation of S17 induced an equilibrium shift toward closed 1, thereby suppressing the binding of p53 to HDM2. This study reveals a multi-state equilibrium of apo-HDM2 and provides new insights into the regulation mechanism of HDM2-p53 interactions.

The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database.

The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.

Simultaneous medium chain fatty acids production and process carbon emissions reduction in a continuous-flow reactor: Re-understanding of carbon flow distribution.

Due to its wide application and high value, the production of medium chain fatty acids (MCFAs) from waste biomass has become one of the worldwide research hotspots. Increasing the carbon element participation from short-chain fatty acids to the form of MCFAs is also conductive to reduce the release of biogas from biological treatment process, because carbon is in the form of MCFAs instead of biogas which directly contribute to process carbon emissions reduction. However, many barriers limiting MCFAs production and application remain to be resolved. Aiming continuous MCFAs production from lactate-rich waste biomass, this study optimized the operation conditions and clarified the main limiting factors and possible mechanisms. The maximum caproic acid concentration of 2.757 g/L were obtained at the Upflow Velocity (ULV) of 1.15 m/h and pH 4.9-5.1. Caproiciproducens, Pseudoramibacter, norank_f_Eubacteriaceae, and Oscillibacter were identified to be the dominant microbial genus responsible for MCFAs production from lactate. The reduction of carbon emissions calculation was also studied in the present processes.

Effect of substrate structure on medium chain fatty acids production and reactor microbiome.

The medium chain fatty acids (MCFAs) produced from organic wastes can replace part fossil-fuel-based products to promote the sustainable development of economy and environment. However, the selection and collocation of feedstocks for MCFAs production are lack of reference basis. This study thereby aimed to investigate how the commonly used electron donor (ED) and substrate configuration affect MCFAs synthesis and then obtain the optimal substrate composition. It was found that the optimized ratios for ethanol/acetate, lactate/acetate, and ethanol/lactate/acetate were 3/1, 2/1, and 2/1/1, respectively, and the optimal substrate concentration was 400 mM C. Combining ethanol and lactate as co-EDs effectively concentrated substrate-carbon-flow (increased by 20-28% than sole ED) on MCFAs synthesis by promoting the elongation of butyrate and reutilization of by-products. As a result, the higher MCFAs yield of 646.22 mg COD/g COD and selectivity of 67.72% were obtained from co-EDs than those from sole ED. Moreover, the key functional bacteria enriched under different ED were also discrepant, which were Clostridium sensu stricto for ethanol, Corynebacterium for lactate, and Veillonella and Oscillibacter for ethanol-lactate, respectively. This study provided a basic but significant reference for the scale-up MCFAs production.

Unraveling the Coupling between Conformational Changes and Ligand Binding in Ribose Binding Protein Using Multiscale Molecular Dynamics and Free-Energy Calculations.

Conformational changes of proteins upon ligand binding are usually explained in terms of several mechanisms including the induced fit, conformational selection, or their mixtures. Due to the slow time scales, conventional molecular dynamics (cMD) simulations based on the atomistic models cannot easily simulate the open-to-closed conformational transition in proteins. In our previous study, we have developed an enhanced sampling scheme (generalized replica exchange with solute tempering selected surface charged residues: gREST_SSCR) for multidomain proteins and applied it to ligand-mediated conformational changes in the G134R mutant of ribose-binding protein (RBPG134R) in solution. The free-energy landscape (FEL) of RBPG134R in the presence of a ribose at the binding site included the open and closed states and two intermediates, open-like and closed-like forms. Only the open and open-like forms existed in the FEL without a ribose. In the current study, the coupling between the conformational changes and ligand binding is further investigated using coarse-grained MD, multiple atomistic cMD, and free-energy calculations. The ribose is easily dissociated from the binding site of wild-type RBP and RBPG134R in the cMD simulations starting from the open and open-like forms. In contrast, it is stable at the binding site in the simulations from the closed and closed-like forms. The free-energy calculations provide the binding affinities of different structures, supporting the results of cMD simulations. Importantly, cMD simulations from the closed-like structures reveal transitions toward the closed one in the presence of a bound ribose. On the basis of the computational results, we propose a molecular mechanism in which conformational selection and induced fit happen in the first and second halves of the open-to-closed transition in RBP, respectively.

The role of gelsolin domain 3 in familial amyloidosis (Finnish type).

In the disease familial amyloidosis, Finnish type (FAF), also known as AGel amyloidosis (AGel), the mechanism by which point mutations in the calcium-regulated actin-severing protein gelsolin lead to furin cleavage is not understood in the intact protein. Here, we provide a structural and biochemical characterization of the FAF variants. X-ray crystallography structures of the FAF mutant gelsolins demonstrate that the mutations do not significantly disrupt the calcium-free conformations of gelsolin. Small-angle X-ray-scattering (SAXS) studies indicate that the FAF calcium-binding site mutants are slower to activate, whereas G167R is as efficient as the wild type. Actin-regulating studies of the gelsolins at the furin cleavage pH (6.5) show that the mutant gelsolins are functional, suggesting that they also adopt relatively normal active conformations. Deletion of gelsolin domains leads to sensitization to furin cleavage, and nanobody-binding protects against furin cleavage. These data indicate instability in the second domain of gelsolin (G2), since loss or gain of G2-stabilizing interactions impacts the efficiency of cleavage by furin. To demonstrate this principle, we engineered non-FAF mutations in G3 that disrupt the G2-G3 interface in the calcium-activated structure. These mutants led to increased furin cleavage. We carried out molecular dynamics (MD) simulations on the FAF and non-FAF mutant G2-G3 fragments of gelsolin. All mutants showed an increase in the distance between the center of masses of the 2 domains (G2 and G3). Since G3 covers the furin cleavage site on G2 in calcium-activated gelsolin, this suggests that destabilization of this interface is a critical step in cleavage.

Exploration of the folding dynamics of human telomeric G-quadruplex with a hybrid atomistic structure-based model.

Structure-based models or Go-like models, which are built from one or multiple particular experimental structures, have been successfully applied to the folding of proteins and RNAs. Recently, a variant termed the hybrid atomistic model advances the description of backbone and side chain interactions of traditional structure-based models, by borrowing the description of local interactions from classical force fields. In this study, we assessed the validity of this model in the folding problem of human telomeric DNA G-quadruplex, where local dihedral terms play important roles. A two-state model was developed and a set of molecular dynamics simulations was conducted to study the folding dynamics of sequence Htel24, which was experimentally validated to adopt two different (3 + 1) hybrid G-quadruplex topologies in K+ solution. Consistent with the experimental observations, the hybrid-1 conformation was found to be more stable and the hybrid-2 conformation was kinetically more favored. The simulations revealed that the hybrid-2 conformation folded in a higher cooperative manner, which may be the reason why it was kinetically more accessible. Moreover, by building a Markov state model, a two-quartet G-quadruplex state and a misfolded state were identified as competing states to complicate the folding process of Htel24. Besides, the simulations also showed that the transition between hybrid-1 and hybrid-2 conformations may proceed an ensemble of hairpin structures. The hybrid atomistic structure-based model reproduced the kinetic partitioning folding dynamics of Htel24 between two different folds, and thus can be used to study the complex folding processes of other G-quadruplex structures.

Metal cofactor modulated folding and target recognition of HIV-1 NCp7.

The HIV-1 nucleocapsid 7 (NCp7) plays crucial roles in multiple stages of HIV-1 life cycle, and its biological functions rely on the binding of zinc ions. Understanding the molecular mechanism of how the zinc ions modulate the conformational dynamics and functions of the NCp7 is essential for the drug development and HIV-1 treatment. In this work, using a structure-based coarse-grained model, we studied the effects of zinc cofactors on the folding and target RNA(SL3) recognition of the NCp7 by molecular dynamics simulations. After reproducing some key properties of the zinc binding and folding of the NCp7 observed in previous experiments, our simulations revealed several interesting features in the metal ion modulated folding and target recognition. Firstly, we showed that the zinc binding makes the folding transition states of the two zinc fingers less structured, which is in line with the Hammond effect observed typically in mutation, temperature or denaturant induced perturbations to protein structure and stability. Secondly, We showed that there exists mutual interplay between the zinc ion binding and NCp7-target recognition. Binding of zinc ions enhances the affinity between the NCp7 and the target RNA, whereas the formation of the NCp7-RNA complex reshapes the intrinsic energy landscape of the NCp7 and increases the stability and zinc affinity of the two zinc fingers. Thirdly, by characterizing the effects of salt concentrations on the target RNA recognition, we showed that the NCp7 achieves optimal balance between the affinity and binding kinetics near the physiologically relevant salt concentrations. In addition, the effects of zinc binding on the inter-domain conformational flexibility and folding cooperativity of the NCp7 were also discussed.

Consequences of Energetic Frustration on the Ligand-Coupled Folding/Dimerization Dynamics of Allosteric Protein S100A12.

Allosteric proteins are featured by energetic degeneracy of two (or more) functionally relevant conformations, therefore their energy landscapes are often locally frustrated. How such frustration affects the protein folding/binding dynamics is not well understood. Here, by using molecular simulations we study the consequences of local frustration in the dimerization dynamics of allosteric proteins based on a homodimer protein S100A12. Despite of the structural symmetry of the two EF-hand motifs in the three-dimensional structures, the S100A12 homodimer shows allosteric behaviors and local frustration only in half of its structural elements, i.e., the C-terminal EF-hand. We showed that such spatially asymmetric location of frustration leads to asymmetric dimerization pathways, in which the dimerization is dominantly initiated by the interchain binding of the minimally frustrated N-terminal EF-hands, achieving optimal balance between the requirements of rapid conformational switching and interchain assembling to the energy landscapes. We also showed that the local frustration, as represented by the double-basin topography of the energy landscape, gives rise to multiple cross-linked dimerization pathways, in which the dimerization is coupled with the allosteric motions of the C-terminal EF-hands. Binding of metal ions tends to reshape the energy landscape and modulate the dimerization pathways. In addition, by employing the frustratometer method, we showed that the highly frustrated residue-pairs in the C-terminal EF-hand are partially unfolded during the conformational transitions of the native homodimer, leading to lowing of free energy barrier. Our results revealed tight interplay between the local frustration of the energy landscape and the dimerization dynamics for allosteric proteins.

A novel missense mutation in the ALPL gene causes dysfunction of the protein.

Hypophosphatasia (HP) is a rare genetic disease caused by mutation in the alkaline phosphatase, liver/bone/kidney (ALPL) gene with highly variable clinical manifestations. Efforts have been made to collect cases with novel mutations and to examine how a missense mutation affects ALPL protein function, which remains difficult to predict. The present study investigated the underlying mechanism of ALPL dysfunction in a patient diagnosed with HP. Bidirectional sequencing of the ALPL gene was conducted in a 5­year­old Chinese girl preliminary diagnosed with childhood HP. Sorting Intolerant from Tolerant (SIFT) and Polymorphism Phenotyping v2 (PolyPhen­2) tools were used to forecast the impact of the mutation on protein function. Site­directed mutagenesis was performed and transfected into cells to verify the role of the specific mutation. Furthermore, the mechanism of the impact was investigated via all­atom molecular dynamics (MD) simulation. The patient demonstrated a compound heterozygote with two missense mutations in the ALPL gene, p.Trp29Arg and p.Ile395Val. Trp29 and Ile395 were determined to be 'tolerable' by SIFT, whereas they were 'possibly damaging' by PolyPhen­2 in terms of conservation. Additionally, HEK293 cells were transfected with plasmids expressing wild type and/or mutated ALPL. Only 4.1% of ALP activity remained when Trp29 was substituted by Arg, whereas 19.1, 33.7, 50.1 and 7.6% ALP activity remained in cells expressing p.Ile395Val, wild type+p.Trp29Arg, wild type+p.Ile395Val and p.Trp29Arg+p.Ile395Val substitutions, respectively. All­atom MD simulation demonstrated that the N­terminal helix of mutated ALPL, where Trp29 is located, separated from the main body of the protein after 30 nsec, and moved freely. These results demonstrated that p.Trp29Arg, as a novel missense mutation in the ALPL gene, reduced the enzymatic activity of ALPL. This effect may be associated with an uncontrolled N­terminal helix. These results provide novel information about the genetic basis of HP, and may facilitate the development of future therapies.

Piecing Together the Allosteric Patterns of Chaperonin GroEL.

Despite considerable efforts, elucidating the allostery of large macromolecular assemblies at a molecular level in solution remains technically challenging due to its structural complexity. Here we have employed an approach combining amide backbone hydrogen/deuterium exchange coupled with mass spectrometry, fluorescence spectroscopy, and molecular simulations to characterize allosteric patterns of chaperonin GroEL, an ∼800 kDa tetradecamer from E. coli. Using available crystal structures of GroEL, we quantitatively map out GroEL allosteric changes in solution by resolving exchange behaviors of 133 overlapping proteolytic peptides with more than 95% sequence coverage. This comprehensive analysis gives a refined resolution down to five residues to pilot the GroEL allosteric determinants, of which the localized dynamics is monitored by tryptophan-mutated GroEL. Furthermore, the GroEL conformational transition is evaluated by molecular dynamics simulations with an atomic-interaction-based coarse-grained model. Collectively, we provide a practical methodology to analyze GroEL allostery in solution.